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Abstract. In the standard treatment of particle oscillations the mass eigenstates are 
implicitly assumed to be scalars and, consequently, the spinorial form of neutrino wave 
functions is not included in the calculations. To analyze this additional effect, we discuss 
the oscillation probability formula obtained by using the Dirac equation as evolution 
equation for the neutrino mass eigenstates. The initial localization of the spinor state 
also implies an interference between positive and negative energy components of mass 
eigenstate wave packets which modifies the standard oscillation probability. 

PACS. 03.65. Pm - 14.60.Pq 

1 Introduction 

Since a long time, particle mixing [1] and oscillation [2,3] continue to stimulate interesting and some- 
times fascinating discussions on the many subtleties of quantum mechanics involved in oscillation 
phenomena. Measurements of various features of the fluxes of atmospheric [4] and solar [5,6] neu- 
trinos have provided, in the last years, evidence for neutrino oscillations and therefore for neutrino 
masses and mixing. In particular, it renewed the interest in understanding the derivation of the flavor 
conversion probability formula and in overcoming the main difficulties hidden in the standard theo- 
retical approaches. In particular, an increasing number of theoretical papers have recently questioned 
the validity of the standard plane wave treatment of oscillations by resorting to intermediate [7-9] 
and external [10, 11] wave packet frameworks. 

The standard plane wave treatment [12,13] is certainly the simplest and probably the most intuitive 
way to introduce the oscillation length and to immediately obtain an expression for the oscillation 
probability. In such a formalism, a plane wave is associated with each mass eigenstate. For the two- 
flavor case, the mass eigenstate phase difference is 



A$ = A{ET-pL) 



(1) 



Thus, an initially pure flavor-eigenstate will be modifled with time and distance. The probability for 
flavor transition is usually expressed in terms of the mixing angle 6 and of the relative phase A<f> by 



P{ua -^Vf3) = sin" [26*] sin" 



(2) 



The Lorentz invariant difference of phase is then conventionally evaluated by setting AT = AL = 
and considering, for ultra-relativistic particles, T L and pi,2 ~ -Ei,2, i-e. 



A(S> = T AE - LAp^ L {AE - Ap) 



At 



■L. 



(3) 
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By using such an approximation, one gets the well-known expression [13 



Uf3'.L) = sin^[2^] sin" 




(4) 



In the plane wave formalism, the most controversial point is certainly represented by the derivation of 
formulas containing extra factors in the oscillation length [14-18]. The use of wave packets allows us to 
understand the origin of these extra factors. In the plane wave approach, it is implicitly assumed that 
at creation the flavor-eigenstate is unique even up to the phase at all points and times of creation. In 
the wave packet treatment, at time T and at a fixed position in the overlapping region, one experiences 
the interference between space points whose separation at creation is given by Aw T and this implies 
that an additional initial phase is automatically included in the wave packet formalism [7,9]. The final 
result contains the difference of phase given in Eq.®. We do not intend here to re-discuss the many 
controversies in the plane wave derivations of the oscillation probability formula. We only remark 
that a plane wave approach leads to conceptual difficulties and fails to explain fundamental aspects 
of particle oscillations (i.e. localization and coherence length). Wave packets eliminates some of these 
problems [19]. In fact, the use of wave packets for propagating mass eigenstatcs {intermediate wave 
packet model) guarantees the existence of a coherence length, avoids the ambiguous approximations in 
the plane wave derivation of the phase difference and, under particular conditions of minimal slippage 
recovers the oscillation probability given in Eq.Q). Unfortunately, it is not easy to determine the 
size of the wave packets at creation and it is not clear whether it makes sense to consider a unique 
time of creation [9,20]. A common argument against the intermediate wave packet formalism is that 
oscillating neutrinos are neither prepared nor observed. Consequently, it would be more convenient 
to write a transition probability between the observable particles involved in the production and 
detection process. This point of view characterizes the so-called external wave packet approach [10, 11] . 
The oscillating particle, described as an internal line of a Feynman diagram by a relativistic mixed 
scalar propagator, propagates between the source and target {external) particles represented by wave 
packets. The function which represents the overlap of the incoming and outgoing wave packets in the 
external wave packet model corresponds to the wave function of the propagating mass eigcnstate in 
the intermediate wave packet formalism. Remarkably, it could be shown that the probability densities 
for ultra-relativistic stable oscillating particles in both frameworks are mathematically equivalent [11]. 
However, the intermediate wave packet picture brings up a problem, as the overlap function takes into 
account not only the properties of the source, but also of the detector. This is unusual for a wave 
packet interpretation and not satisfying for causality [11]. This point was clarified by Giunti [10] who 
solves this problem by proposing an improved version of the intermediate wave packet model where 
the wave packet of the oscillating particle is explicitly computed with field-theoretical methods in 
terms of external wave packets. Despite not being applied in a completely free way, the {intermediate) 
wave packet treatment commonly simplifies the discussion of some physical aspects going with the 
oscillation phenomena [9,17]. Thus, it makes sense, as a preliminary investigation, to consider a wave 
packet associated with the propagating particle. 

In this paper, we aim to investigate how the oscillation formula is modified by using fermionic 
instead of scalar particles. To do it, we shall use the Dirac equation as the evolution equation for 
the mass eigenstates. Before introducing the Dirac formalism, in section II, we briefly review the 
intermediate wave packet model for scalar particles [19]. In this section, by choosing a gaussian wave 
packet to describe the localization of our initial flavor state, we obtain an analytical expression for the 
flavor conversion probability. This allows us to identify the wave packet slippage and spreading effects. 
In section III, we introduce the Dirac formalism and show that a superposition of both positive and 
negative frequency solutions of Dirac equation is often a necessary condition to correctly describe 
the time evolution of the mass eigenstate wave packets. We give, for strictly peaked momentum 
distributions and ultra-relativistic particles, an analytic expression for the Dirac flavor conversion 
probability. The results obtained in the context of a wave packet treatment of oscillation phenomena 
are (briefly) compared with quantum field theory calculations [11,21,22]. This allow to understand how 
our analysis could be included within the external wave packets formalism. We draw our conclusions 
in Section IV. 
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2 Gaussian wave packets 

The main aspects of oscillation phenomena can be understood by studying the two flavor problem. In 
addition, substantial mathematical simplification results from the assumption that the space depen- 
dence of wave functions is one-dimcnsional (z-axis). Therefore, we shall use these simplifications to 
calculate the oscillation probabilities. In this context, the time evolution of fiavor wave packets can 
be described by 

<I>(z,t) = (f>i{z,t) cos 9 v>i + 02(-z, sinfl 1/2 
— [(j)i(z,t) cos^ 6 + 4)2{z,t) sin^ 6] 

+ [01 (z, t) — (/)2 (z, t)] cos Q sin Q vp 
= (j)a{z,t;9)i>a + <l>i3{z,t;e)v>p, (5) 

where i^a and are flavor-eigcnstatcs and and fa ai'c mass eigenstates. The probability of finding 
a flavor state at the instant t is equal to the integrand squared modulus of the vp coefficient 

P,c.,..K 1^/3; t) = J^^ dz \^f3{z,t;0)\' 

= 2i4lM{l-lNT_,„.(t)}, 



(6) 



where lNT^^^i^^{t) represents the interference oscillating term between the (scalar) mass eigenstate 
wave packets 0i(z,t) and 02 (-2,^), i-c 



/_ dz(j)\{z,t)(l)2{z,t) 



lNT,,„,„,(t) = Re 

Let us consider mass eigenstate wave packets given at time i = by 

1 _ _ 

2 



exp 



exp [ipi 



The wave functions which describe their time evolution are 

„+ 



where 
and 



(j)i{z,t) = ^(p{p2.-Pi)exp[-~-iE.i{p^)t + ipzz], 



(7) 



(8) 



(9) 



viPz -Pi) = (2™^)-* exp 



{Pz -PiYa" 



In order to obtain the oscillation probability, wc must calculate the interference term lNT3c„,„r(t), i. 
e. we have to solve the following integral 



/_ ^ fiPz - Pi)v{Pz - P2) exp [-Z AE{pz) t] = 



exp 



JaApy 



L ^v'iPz~p)cxp[~^AEipz)t], 



(10) 



where we have changed the z integration into a pz integration and introduced the quantities Ap = 
Pi — Pi, P = ^{Pi + P2) and AE{pz) = E^^p^) — E^iPz)- The oscillation term is bounded by the 
exponential function of a Ap at any instant of time. Under this condition wc could never observe a pure 
fiavor-eigenstate. Besides, oscillations are considerably suppressed if a Ap > 1. A necessary condition 
to observe oscillations is that aAp ^ 1. This constraint can also be expressed by 5p ^ Ap where 
5p is the momentum uncertainty of the particle. The overlap between the momentum distributions is 
indeed relevant only for 5p 3> Ap. Consequently, without loss of generality, we can assume 



Int^ 



Re 



L ^^HPz-p)exp[-tAE,{pz)t] 



(11) 



4 



Alex E. Bernardini, Stefano De Leo: Dirac spinors and flavor oscillations 



In litterature, this equation is often obtained by assuming two mass eigenstate wave packets described 
by the "same" momentum distribution centered around the average momentum p{= p„). This sim- 
phfying hypothesis also guarantees instantaneous creation of a pure flavor eigenstate at t = 0. In 
fact, for (j}i{z,0) = 4>2{z,0) we get from Eq.© 



cxp [ip z] 



and 

To analytically solve the integral in Ea. lfTT)) . let us rewrite the energy Ei[pz) as follows 



(12) 
(13) 



E^{p.)=P 

= p(i + x) 



P 



P 



1 + 



where 



(1+x) 
and 



X 



p 



(14) 
(15) 



In what follows, we shall consider ultra-relativistic particles and assume a sharply peaked momentum 
distribution, i. e. 

mi<p => Ci < 1 and dp ^ p X < 1- 

Let us now expand the energy Ei{pz) in a power series of Q and We choose to cut off the power 
series terms of order Q 



so that 

E^{p^) 



P 



= P 

In this case, the energy difference becomes 

Ai?(p.; 



1 + 



2(1+X) 



oo 



(16) 



(17) 



By considering only the first term in the x expansion, we reproduce the plane wave result. Indeed, 

Ai?["l(p,)=p^. (18) 

An approximation of order ^ 1) in Eq. ljlTI) requires some constraints on x- Since we cut off 

terms of order A^^ and we wish to consider terms up to x'^A^' in Ea. H17() . it is necessary to satisfy the 
constraint x^^C > '^hich implies x > (C = ^^"5^)- same time, for eliminating x'^'''^AC, 

we have to impose x'^~'^^AC < which can be rewritten as x < C^^- In this way, an approximation 

of order x*^ will be consistent in the range < x ^ . Meanwhile, the integral in Ea. Hll|) can 
be solved analytically only when k < 2. By taking into account terms up to the order x^, the energy 
difference becomes 

Ai?Pl(p,)-p-^(l-X + x'). (19) 
If we substitute (|19|l in Ea. (|ll|) we obtain 

lNT,,<.,_(t) « 

^^{7tl-2dx cxp 



X exp 



IS{t) + ^Q{t) 



-m)]} 
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where 



S{t) = ^t, Q{t) = -jS^ and 7^(^) 



Am^t 



(21) 



By suppressing the variable (i) dependence, we can rewrite lNT^^„,„r(t) as 

lNT„„,„^(i) w exp 



'4(l+K^) 



(1 + ^2)7+1 

2(l+7?,2) 



4(l+7?,2) 



(1+^2)^-1 ■ 



(22) 



The above result deserve some comments. The wave packet spreading is parameterized by TZ(t). At 
the same time, the slippage effect between mass eigenstates is predominantly quantified by the Q{t) 
parameter. The spreading of wave packets is a secondary effect with respect to the slippage since from 
Ea. (|21|l we can write 

Sg - ^- (23) 
Under minimal spreading conditions, i. e. when TZ <^ 1, the Ea. (|22() becomes 



Int 



exp 



cos 



g^(i-R^) 

4 



TZ ■ 

sm 



}■ 



(24) 



where the oscillating character is predominantly given by the cosine function behavior. The exponential 
term with TZ(t) extends the interference between the mass eigenstatc wave packets, and consequently 
the oscillating character, for (a little) longer times. Taking into account terms up to the order x in 
the Ea. l(T7|l . we can write 



Ai?W(p,)=p^(l-x) 
and compute the oscillation probability with the leading corrections due to the slippage effect, 



■■[20] 



1 — exp 



Asm], 



(25) 



(26) 



which corresponds to the same result obtained by [9]. Under minimal slippage conditions, i. e. when 
Q{t) <^ 1, the Ea. H26|) reproduces the plane wave formula 



sin'' [26)] 
2 



cos 



[St]] 



2iHlIM{i_cos [S{t)]} 



sin' \2e] sin' 



4p 



(27) 



3 Dirac formalism 

The results in the previous section have been obtained by considering scalar mass eigenstates. Neutri- 
nos are, however, fermions. The time evolution of a spin one-half particle have to be described by the 
Dirac equation. To introduce the fermionic character in the study of quantum oscillation phenomena, 
we shall use the Dirac equation as the evolution equation for the mass eigenstates. The Eq.(|Sll now 
becomes 

\E'(z, t) = ipi{z, t) cos 6* i/i + ip2{z, t) sin 6 v-^ 
= [tlJi{z, t) cos' 9 + ip2{z, t) sin' 9] 

+ [i^i[z, t) — ip2{z, t)] cos 6* sin 6* i>'/3 
= ^^[z,t-9)uo,+'4>fi{z,t-9)up, (28) 
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where ipi{z,t) satisfies the Dirac equation for a mass mi. The natural extension of Ea. H13|l reads 

Mz,0,e)=cj)^{z,0,9)w (29) 
where w is a constant spinor which satisfies the normalization condition w^w = 1. 

3.1 Dirac wave packets and the oscillation formula 

To describe the time evolution of mass eigenstatc Dirac wave packets, we could be inclined to superpose 
only positive frequency solutions of the Dirac equation. It seems, at first glance, a reasonable choice. 
However, when the initial state has the form given in Ea. H29|) . it is necessary to superpose both 
positive and negative frequency solutions of Dirac equation. Let us clear up this point. The flavor 
state ipaiz, t, 9) is now expressed in terms of 

ip^{z, t) = ^ cxp [ipzz] 

+ dt*{-PzXi-Pz) exp [+tE,{pz)t]}. (30) 

At time t = the mass eigenstatc wave functions satisfy ■01(^1 0) ~ ip2(z,0) (this guarantees the 
instantaneous creation of a pure flavor-cigenstate Va we have appointed in section II) . The Fourier 
transform of ipi^z, 0) is 

[Kip.MiPz) + dn-p.xi-p.)] . (31) 

s=l,2 

By observing that the Fourier transform of (/)a{z, 0, 9) is given by (p{pz — p) (see Ea. (|13|l ). we imme- 
diately obtain the Fourier transform of ij^a^z, 0, 9), 

^{p. -p)w=Y IMiPzKiP.) + dn-pMi-Pz)] ■ (32) 

s=l,2 

Using the orthogonality properties of Dirac spinors, we find [23] 

KijPz) = ^p{Pz -p)uf{pz)w, 

dTi-Pz) = v{Pz~p)vf{-p,)w. (33) 

These coefficients carry an important physical information. For any initial state which has the form 
given in Ea. (|29|l . the the negative frequency solution coefficients df{~pz) necessarily provides a non- 
null contribution to the time evolving wave packet. This obliges us to take the complete set of Dirac 
equation solutions to construct the wave packet. Only if we consider a momentum distribution given 
by a delta function (plane wave limit) and suppose an initial spinor w being a positive energy mass 
eigenstatc with momentum p, the contribution due to d'l*{-pz) will be null. 

Having introduced the Dirac wave packet prescription, we are now in a position to calculate the 
flavor conversion formula. The following calculations do not depend on the gamma matrix represen- 
tation. By substituting the coefficients given by Ea. (|33|l in Ea.H30|l and using the well-known spinor 
properties [23], 

UdPz)u,[p^) = — 



s=l,2 



SI ^-s, N l°Ei{p^)+-i^p^-m, 
Y V, i-PzK i-Pz) = , (34) 
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we obtain 

dp 

-^viPz- P) exp [ipzz]{cos [E.i{pz)t] 

»7° il'Pz + mi) . 
sin[£;,(p^)i]}w. (35) 

By simple mathematical manipulations, the new interference oscillating term will be written as 

lNT,_(t)=/;~^^^(p,_p) 

X {(l-F(p,))cos [AEip,)t] 

+ Fip,)cos[2E{p,)t]} (36) 



where 
and 



Erjp,) + E,{p,) 

2 



1 pI+ m^m^ 
F(Pz) = 77 - 



2 2E,{p,)E,{p,) 



What is interesting about the result in Ea. (|36|l is that it was obtained without any assumption on the 
initial spinor w. Otherwise, the initial spinor carries some fundamental physical information about 
the created state. And this could be relevant in the study of chiral oscillations [28] where the initial 
state plays a fundamental role. With respect to the standard treatment of neutrino oscillations done 
by using scalar wave packets and leading to the interference term 

INT..„,„(0 = / ^ip^Pz- p) COS [\E{pz) <], (37) 

we note in l¥iToir-aa{i) two additional terms. In the first one, the stonrfani oscillating term cos [A_B(pz) t], 
which arises from the interference between mass eigenstate components of equal sign frequencies, is 
multiplied by a new factor obtained by the products 

u\{pz)u2{Pz), vli-Pz)v2i-Pz) andh.c. 

The second one, cos[2E{pz)t\, is a new oscillating term which comes from the interference between 
mass eigenstate components of positive and negative frequencies. The factor multiplying such an 
additional oscillating term is obtained by the products 

uliPz)v2i-Pz), v\{-pz)u2{pz) andh.c. 

The new oscillations have very high frequencies. Such a peculiar oscillating behavior is similar to 
the phenomenon referred to as Zitterhewegung. In atomic physics, the electron exhibits this violent 
quantum fluctuation in the position and becomes sensitive to an effective potential which explains 
the Darwin term in the hydrogen atom [24]. We shall see later that, at the instant of creation, such 
rapid oscillations introduce a small modification in the oscillation formula. 

We plot the function F{pz) in Figure ^ We can readily observe that it goes rapidly to zero for 
Pz 3> TTii 2, it has a minimum sX pz = Q and two maxima at pz = iy^mTrol. The maximum value of 
F{pz) is ' 

Fma.{pz) = -[l- ^ : - ] 38 

2 \ mi + m.2 J 

which vanishes in the limit ~ m^. As we can see in Figure ^ the new effects are relevant only 
when Am mj ^ mj. 
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Pi ^ '"i 

Fig. 1. The function F{p^) is plotted for different values of the ratio between mi and m^. For a momentum 
distribution sharply peaked around p S> "ii,2, F{pz) does not play a significant role in the "modified" oscil- 
lation formula. In the case of mi ~ m2, independently of the value of p and of the momentum distribution 
width, the maximum values of F{pz) are negligible and consequently F{pz) is practically suppressed in the 
calculation (see amplification in the upper box). 



3.2 The oscillation formula with spreading 

To quantify the new effects exhibited in the oscillation probability formula let us calculate the integral 
of Eq. (|36|l by assuming an ultra-relativistic particle and following the localization condition given by 
(|16|l . We use the same criteria adopted in section II to expand the energy H14() in a power series of d 
and X- A second order approximation in x allows us to write the energy dependent terms as 



(39) 



2pi J 2p~X 



(40) 

(4f)'(1-2x + 3x^). (41) 

Since we have approximated not only AE{pz), but also F{pz) and E{pz), the range of validity for an 

analytical approximation of order x'' is now given by ^■^^ < X 1^ C'h^. By substituting ip{pz — p) 
and the approximations 1)39141(1 in Eg. ((36(1 . we obtain 



lNTo.,„(t) 



ap 



(27r)5 



r /_ dx exp 



{[i-(4f)'(1-2x + 3x^) 

^"'*(i-x + x^)' 



X cos 



2p 



+ (lf)'(l-2x + 3x^) 
X cos [pt (2(1 + x) + - X + X'))] } 

A new integrating variable a = and the coefficients 



(42) 



(43) 
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enable us to write Ea. H42(l in the form 

lNT^™(t) « ^ exp [-a^] 



6 



X exp[-i5_(t) -^Q_(^)cr-^7^_(^)cr'] 



X exp - iQ+it)<T - ^7^+(^)o•']} 

= Re [H_{t)G_ [t) + H^{t)G+{t)] 



where 



G±(t) 



i+in±(t) 

are obtained in the same way as (|26(l . and 



exp 



g±(«) 

4(l+^7^±(t)) 



^ 1 All 



arise from the new coefficients which include F{pz). 



(44) 
(45) 



(46) 



3.3 The oscillation formula without spreading 

A more satisfactory interpretation of the modifications introduced by the Dirac formalism is given 
when we restrict our study to a first order approximation in x, i. e. without considering the wave 
packet spreading. In fact, we could take into account terms up to the order x in the Eas. H39l4f |l and 
obtain a simpler approximation, 

l^T^^Ut) « C ^ { [l - (t^)' (l - ^-)^ 
X exp [-iS_ (t) - iQ_ (t)cr - cr"] 

X exp [~iS+{t) - iQ+{t)a - a']} 



X exp 



QUt) 

4 

ap 



X exp 

By using the explicit expressions for Q±{t) and S±{t) we get 



exp 



2p 



(47) 
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As we have already noticed, the osciUating functions going with the second exponential function in 
Ea. (|48|l arise from the interference between positive and negative frequency solutions of the Dirac 
equation. It produces very high frequency oscillations which is similar to the quoted phenomenon of 
Zitterbewegung [24]. The oscillation length which characterizes the very high frequency oscillations is 
given by L^^^ ~ Obviously, L^^^ is much smaller than the standard oscillation length given 

by sc = Sn^ • means that the propagating particle exhibits a violent quantum fluctuation of its 
flavor quantum number around a flavor average value which oscillates with iosc ■ Meanwhile, except 
at times t 0, it provides a practically null contribution to the oscillation probability. To explain 
such a statement, let us suppose that an experimental measurement takes place after a time T w L 
for ultra-relativistic particles. The observability conditions impose that the propagation distance L 
must be larger than the wave packet localization a. Since the (second) exponential function vanishes 
when L ^ a, for measurable distances, the effective flavor conversion formula will not contain such 
very high frequency oscillation terms, and can be written as 



X |l — exp 



Ar: 



2^/2ap^ 



COS 



Am 
2p 



2p 



2p 

}}■ 



2-i 



(49) 



For distances which are restrict to the interval a L <^ we observe the minimal slippage 

between the wave packets. In this case, we could suddenly approximate the oscillation probability to 



sin^ [26] 
2 



Y \2V2ap^) 
1 - 



Am 
2p 



cos 



(50) 



however, we reemphasizc that it is not valid for T ^ L ^ when the rapid oscillations are still relevant 
{L < a). By comparing the result of Eq. (|5n|l with the scalar oscillation probability of Eq. , we notice 

a deviation of the order (^-^^^ that appears as an additional coefficient of the cosine function. It is 

not relevant in the ultra-relativistic limit as we have noticed after studying the function F{pz). 



3.4 A brief extension to quantum field treatment 

To finalize our study, we try to establish a tenuous correspondence between our results and the QFT 
treatment. It was extensively demonstrated in the literature [10, 20, 21] that the oscillating particle 
cannot be treated in isolation. The oscillation process must be considered globally: the oscillating states 
become intermediate states, not directly observed, which propagate between a source and a detector. 
This idea can be implemented in QFT when the intermediate oscillating states are represented by 
internal lines of Fcynman diagrams and the interacting particles at source/detector are described 
by external wave packets [11,21]. In this context, let us consider the weak flavor-changing processes 
occurring through the intermediate propagation of a neutrino. 

Pi —t Pp + a + i^a (oscillation) vp + Dj f3 + Dp (51) 

where Pi and Pp (Dj and Dp) are respectively the initial and final production (detection) particles. 
The amplitude for the process is represented by 

{Pp,Dp\T {exp[-i J dx^Tii]) -l\Pi,Di) (52) 

where Hj is the interaction Hamiltonian for the intermediate particle and T is the time ordering 
operator. After some mathematical manipulations [11], this amplitude can be represented by the 
integral 

X G{E, p,tD,tp) exp[ip ■ (xd ~ xp)] (53) 
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where the function F{E,p) represents the overlap of the incommg and outgoing wave packets, both 
at the source and at the detector, and the Green function in the momentum space, G{E,p,t]j,tp), 
represents the fermion propagator which carries the information of the osciUation process. The overlap 
function is independent of production and detection times and positions {tp, to, xp, x^i) and depends 
on the the directions of incoming and outgoing momenta. In certain way, the physical conditions of 
source and detector, in terms of time and space intervals, are better defined in this framework than 
in the intermediate wave packet framework. Anyway, to understand the oscillation process we must 
turn back to the definition of mixing in quantum mechanics. It is similar in field theory, except that 
it applies to fields, not to physical states. This difference allows to bypass the problems arising in 
the definition of flavor and mass bases [11]. In one-dimcnsional spatial coordinates, the mixing is 
illustrated by the unitary transformation 

^l^,{z,t;9) ^ g-\e-t)ij,{z,t)g{e;t) (54) 

as the result of the noncoincidence of the flavor basis (cr = a, (3) and the mass basis (i = 1, 2). The 
Ea. (|54|l gives the the Ea. (|28|l when the generator of mixing transformations g{9\ t) is given by 

^(e';t) = exp[6' j dz^^{z,t)'il;^{z,t) 

(55) 

By taking the one-dimensional representation of Ea. H53|l . the propagator G{E,p2,ti:),tp) can also be 
written in the flavor basis as 

G''"{e; E,p,,T) = g-\e- 1) g{e,p,,t) g{e- 1) 

= g-\e; t) G{E,p,,tD,tp) g{9; t) (56) 

with T ^ tp) — tp. 

In particular, by following the Blasone and Vitiello (BV) prescription [22,25], the definition of a 
Fock space of weak eigenstates becomes possible and a nonperturbative fiavor oscillation amplitude 
can be derived. In this case, the complete Lagrangian (density) is split in a propagation Lagrangian, 

£p(z, t) = Vii(z, t) [i^ - TOi) V'i(2, t) 

+ i},iz,t) {ip ~m,) Mz,t), (57) 

and an interaction Lagrangian 

Ci{z, t) ^ ipaiz, t; 9) {ip ~ ma) -ipa^z, t; 0) 

+ tp0{z,t;0) {ilj) ~mp) ipp{z,t]9) 
-map {'ipa{z,t;6)il;p{z,t;0) 

+ i;p{z,t;0)i;aiz,t]0)), (58) 

where 

i^a{p) = '^1(2) cos^ 6 + 7712(1) sin^ 

and 

fna/i = (ttt-i — TO2) COS sin 0. 

In general, the two subsets of the Lagrangian can be distinguished if there is a flavor transformation 
which is a symmetry of Li{z, t) but not oi Lp{z, t). Particle mixing occurs if the propagator built from 
Cp{z,t), and representing the creation of a particle of flavor a at point z and the annihilation of a 
particle of flavor (3 at point z', is not diagonal, i.e. not zero for (3 = a. The free fields il)i{z,t) can be 
quantized in the usual way by rewriting the momentum distributions bf(jiz) and df*{-pz) in Ea. (|30|l 
as creation and annihilation operators B^{pz) and D'p{-pz). The interacting fields arc then given by 

Mz,t) ^ /^""^ i|^exp[zp,z] ^ {S^(p,;t)<(p,;t) 

s=l,2 

+ D:;i-Pz;t)vU-Pz;t)} (59) 
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where the new flavor creation and annihilation operators which satisfy canonical anticommutation 
relations are defined by means of Bogoliubov transformations [25] as 



and 



Dl{-p,;t) = g-\e-t) DU-Pz)g{0;t) 



By following the BV prescription [22], which takes into account the above definitions, it was demon- 
strated [27] that the flavor conversion formula can be written as 



+ \{D^p{-p;t),B:,{p;t)}\ 



(60) 



which is calculated without considering the localization conditions imposed by wave packets, i. e. by 
assuming that pz w p. When the explicit form of the flavor annihilation and creation operators are 
substituted in Eg. 1)60(1 . it was also demonstrated [25] that the flavor oscillation formula becomes 

2iHlIM{(i„F(p))cos[Ai?(p)i] 

+ F{p) cos [2E{p)t]} 



in^ [20]{[l-(^)^ 



Am 
2p 



ptil 



Am" 
4p 



(61) 



where the last approximation takes place in the rclativistic limit p 3> ^m^m^. After some simple 
mathematical manipulations, the Eg. 1)61(1 gives exactly the oscillation probability PowaX^oL i^p', L) 
calculated from Eg. ((48(1 when it is assumed that the wave packet width a tends to infinity and t ^ L. 

This new oscillation formula tends to the standard one 10} in the ultra-relativistic limit. If the 
mass eigenstates were nearly degenerate, we could have focused on the case of a nonrelativistic os- 
cillating particle having very distinct mass eigenstates. Under these conditions, the quantum theory 
of measurement says that interference vanishes. Therefore, as we have already appointed, the effects 
are, under realistic conditions, far from observable. Besides, in spite of working on a QFT framework, 
the lack of observability conditions must be overcome by implementing external wave packets, i. e. 
by calculating the explicit form of Eg. 1(53(1 for fermions. Such a procedure was applied by Beuthe for 
scalar particles [11] and, in a very particular analysis, with basis on the BV calculations and on our 
intermediate wave packet results, it could be extended to the fcrmionic case. 



4 Conclusions 

In this paper we have computed the modifications to the flavor conversion probability caused by the 
introduction of the spinorial form of neutrino wave functions. To describe the time evolution of the 
mass eigenstates, we have introduced wave packets constructed by superposing the Dirac equation 
solutions. By following an analytical study with gaussian wave packets we have computed the new 
effects that can be observed in the flavor conversion probability formula. Our study leads to the 
conclusion that the fermionic nature of the particles and the interference between positive and negative 
frequency components of mass eigenstate wave packets modify the standard oscillation probability 
which is obtained by implicitly assuming a scalar nature of the mass eigenstates. Nevertheless, under 
particular assumptions, i.e. ultra-relativistic particles and sharply peaked momentum distributions, 

these modifications introduce correction factors proportional to (^^^^ which are negligible in the 

ultra-relativistic limit. 

We know, however, that the most rigorous treatment of oscillations might be done in the quantum 
field theory framework. Meanwhile, the prescription of oscillating neutrinos as Dirac spinors was 
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not yet completely and accurately described in a quantum field formalism. The BV model [22,26] to 
neutrino/particle mixing and oscillations is the most preeminent trying to this aim. They have attempt 
to define a Fock space of weak eigenstates and to derive a nonperturbativc oscillation formula. Flavor 
creation and annihilation operators, satisfying canonical (anti)comutation relations, arc defined by 
means of Bogoliubov transformations. As a result, new oscillation formulas are obtained for fermions 
and bosons, with the oscillation frequency depending not only on the difference but also on the sum 
of the energies of the different mass eigenstates. 

By using Dirac wave packets, we have reproduced an oscillation probability formula with the same 
mathematical structure as those obtained in the BV model [22, 26] in a QFT framework. The study 
with Dirac wave packets enables us to quantify separately each new effect present in the oscillation 
formula. Imposing the initial constraint where we have a pure flavor-eigenstate at time of creation 
t = for any constant spinor w, we could calculate the contribution of new effects to the oscillation 
probability. Particularly, we have noticed that a term of very high oscillation frequency depending 
on the sum of energies introduces a very small modification in the characteristic of the oscillation 
phenomena. In addition, the spinorial form of the wave functions subtly modifies the coefficients of 
the oscillating terms in the flavor conversion formula. 

To conclude, we emphasize one more conceptual aspect arising from the Dirac formalism. Dirac 
wave packets enable us to develop a study of chiral oscillations [28] . In the standard model of flavor- 
changing interactions, neutrinos with positive chirality are decoupled from the neutrino absorbing 
charged weak currents. In the ultra-relativistic limit, a state with left-handed helicity is practically 
a state with negative chirality. If the interactions at the source and detector are chiral, only the 
component with negative chirality contributes to the propagation. Therefore, the possibility of chiral 
oscillations can subtly modify the oscillation formula. In this context, the study of chiral and flavor 
oscillations could also deserve some further specific studies. 
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